On extreme field limits in high power laser matter interactions: 
radiation dominant regimes in high intensity electromagnetic wave 
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Abstract We discuss the key important 
regimes of electromagnetic field interaction with 
charged particles. Main attention is paid to the 
nonlinear Thomson/ Compton scattering regime 
with the radiation friction and quantum elec- 
trodynamics effects taken into account. This 
process opens a channel of high efficiency elec- 
tromagnetic energy conversion into hard elec- 
tromagnetic radiation in the form of ultra short 
high power gamma ray flashes. 



I. INTRODUCTION 

Radiation of present day lasers approaches 
the intensity regimes where in the electromag- 
netic (EM) wave interaction with matter the 
radiation friction effects on the charged parti- 
cle dynamics become dominant [TJ [2] . At these 
limits the electron dynamics become dissipa- 
tive [3]-[6] with fast conversion of the EM wave 
energy to hard EM radiation, which for typi- 
cal laser parameters is in the gamma-ray range 
[71 [8j . For laser radiation with 1 jam. wavelength 
the radiation friction force changes the scenario 
of the EM wave interaction with matter at the 
intensity of about Ir ^ lO^^W/cm^. For the 
laser intensity close to Ir also novel physics of 
abundant electron-positron pair creation comes 



into play [9l [TO] (see [H [12] and [HITS]). In 

this regime, the electron (positron) interaction 
with the EM field is principally determined by a 
counterplay between the radiation friction and 
quantum effects. The quantum electrodynamics 
(QED) effects weaken the EM emission by the 
relativistic electron resulting in the lowering of 
the radiation friction [17, 18j. In an extremely 
high intensity EM field vacuum looses its prop- 
erty to be a empty substance showng such vari- 
ous nonlinear QED processes as vacuum polar- 
ization, electron-positron pair plasma creation, 
and other properties depending on the EM field 
amplitude. 

In Fig. [l]we present a schematic of nonlinear 
QED processes which realization depends not 
only on the EM radiation intensity, but also on 
the photon and charged particle energy and on 
the EM field configuration [TOHTO]. Fig. [l] il- 
lustrates a transition from the relativistic in- 
teraction regime (Rel) to dominant radiation 
friction (RF) and nonlinear Thomson scatter- 
ing (NTS) through the limit when the quantum 
electrodynamics comes into play (QED) and 
to electron-positron avalance/cascade develop- 
ment (A/G), towards nonlinear vacuum with 
electron-positron creation (Es) and nonlinear 
vacuum polarization (VP) while the EM field 
intensity grows. 
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FIG. 1. Schematic of nonlinear QED processes. 
Laser-matter interaction transition from the rel- 
ativistic regime (Rel) to dominant radiation fric- 
tion (RF) and nonhnear Thomson scattering (NTS) 
through the hmit when quantum electrodynamics 
comes into play (QED) and to electron-positron 
avalance/cascade development (A/C), towards non- 
linear vacuum with electron-positron creation (Es) 
and nonlinear vacuum polarization (VP) while the 
EM field intensity increases. 



below a summation over repeating indices ja = 
0, 1,2,3 is assumed. 

A relativistic electron interacting with an EM 
wave emits high energy photons. Here and be- 
low for the sake of simplicity we analyse the 
dynamics of a radiating electron in a homoge- 
neous rotating electric field, which corresponds 
to the nodes of two colliding EM waves, where 
the wave magnetic field vanishes, and/or to the 
electron interaction with an EM wave in near- 
critical plasmas in the frame of reference moving 
with the group velocity of the wave [19 . 

In the regime of Nonlinear Thomson Scatter- 
ing (NTS) the power emitted is proportional to 
the fourth power of its energy, 7TieC^7, [20] 



' eradrrieC^ujj^. 
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The probabilities of the processes involving 
extremely high intensity EM field interaction 
with electrons, positrons and photons are de- 
termined by several dimensionless parameters. 

When the normalized dimensionless EM wave 
amplitude 



eEX 



(1) 



exceeds unity, a > 1, the energy of the elec- 
tron quivering in the field of the wave becomes 
relativistic. Here A = 27tc/uj with uj being the 
EM wave frequency. The EM wave intensity is 
expressed via the normalized amplitude, a, as 



4^e2 

1.37 X 10-a^ f 1/^ V ^ 
\ X J cm^ 
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In a plane EM wave the parameter a is related 
to the Lorentz invariant, which being expressed 
via the 4-potential of the electromagnetic field. 



proportional to the ratio of the classical electron 
radius rg = e^ /rUeC^ = 2.8 x 10~^^ cm and the 
EM wave wavelength A characterizes the role of 
radiation losses. The maximal rate at which an 
electron can acquire the energy from the EM 
field is approximately equal to rrieC^uja. The 
condition of the balance between the acquired 
and lost energy for the electron Lorentz factor 
equal to 7e = a shows that the radiation effects 



become dominant at a > arad 
the EM wave intensity above 
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The characteristic frequency of the emitted 
radiation is proportional to the cube of the elec- 
tron energy, ujm ~ 0.3co'7g. QED effects be- 
come important, when the energy of the pho- 
ton generated by Thomson (Compton) scatter- 
ing is of the order of the electron energy, i.e. 
hum ^ rrieC^ie' If 7e = <^ this yields the QED 



A^, is equal to a = e^A^A^/rrieC^. Here and limit on the EM field amplitude, o? /as > 1. 



Here the dimensionless parameter 



as 



eEsX 



2^ 



27rmeC^ huj 
4.2 X 10^ 



1/im 
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is the normaUzed critical electric field of quan- 
tum electrodynamics, Es = mlc^/eh [21 , with 
Ac = 27rh/meC = 2.42 x 10"^°cm being the 
Compton wavelength. The EM radiation inten- 
sity for the wave with the amplitude as is 



Is 



Ane^h? 



2.36 xlO'^^. (7) 
cm^ 



The above obtained QED limit, a^ /as > 1, cor- 
responds to the condition Xe > 1, where rela- 
tivistic and gauge invariant parameter Xe 



Xe 



EsTUeC 



a 

as 



(8) 



characterizes the probability of the gamma- 
photon emission by the electron with 4- 
momentum p^, in the field of the EM wave. The 
4-tensor of the EM field is defined as F, 



jlV 



d^Ay — dyA^. The QED limit is reached for the 
EM wave intensity of the order of 
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The dimensionless parameter 



X7 



EsTUeC 



tto U 



(10) 



determines the probability of the electron- 
positron pair creation by the photon with the 
energy huj^ in the EM field via the Breit- 
Wheeler process [22] |23] . 

Nonlinear vacuum properties are determined 
by the Poincare invariants. 



and 



F = (B^ - e2)/2 



(E-B). 
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FIG. 2. Regimes of EM field interaction in 

the space of parameters E/^^, x^, and Xe- The 
limit K/Es -^ 1 corresponds to the nonlinear QED 
vacuum regimes with the electron-positron pair 
ceation from vacuum and photon-photon interac- 
tion. When the parameter Xe becomes large, multi- 
photon Compton scattering results in the high en- 
ergy photon generation. For X7 > 1 the mutiphoton 
Breit- Wheeler process results in electron-positron 
pair generation via the gamma-photon interaction 
with a strong EM field. The nonlinear Thomson 
scattering regime is realized for a ^ 1 with the ra- 
diation friction effects coming in play at a > arad- 



The probability of the electron-positron pair 
creation depends on 



E 



Vf2 + G2 _ if 



and 



VF2 + G2 + F, 



(13) 



(14) 



which are the electric and magnetic fields in the 
frame of reference where they are parallel. 

Figure [2] illustrates different regimes of the 
EM interaction in the space of parameters 
¥,/Es^ X77 ^^d Xe- In the limit 'E/Es -^ 1 
the EM waves can create electron-positron pairs 
from vacuum [2, 25 , the EM wave can interact 
via photon-photon collisions [26] and vacuum 
polarization [27 . When the parameter Xe be- 
comes large the multiphoton Compton scatter- 
ing results in the high energy photon genera- 
tion. For X7 > 1 the mutiphoton Breit- Wheeler 



process results in electron-positron pair gener- 
ation via the gamma-photon interaction with 
the strong EM field. The nonlinear Thomson 
scattering (NTS) regime is realized for a ^ 1 
with the scattering cross section depending on 

— 1/3 

the EM field amplitude. If a > arad = ^rad 
radiation friction effects play a key role. 

The comparision of expressions given by Eqs. 
(pi and Q shows that for the EM wave length 
equal to A ?^ 0.8/im the intensities Ir and Iq 
are of the same order of magnitude as has been 
noted in Ref. [28 . The curves Ir{oj) and Iq{oj) 
intersect to each other at the frequency equal to 



UJi 



4 

18^3 



(15) 



corresponding to the wavelength of Ai = 
0.821/im and the photon energy of the order of 
1.5 eV. 

It is convenient to rewrite expressions ([5| and 
(|9| for Ir{uj) and Iq{oj) in terms of the fre- 
quency cji as 
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In the present paper we mainly pay attention 
to the nonlinear Thomson/Compton scattering 
regime when the radiation friction and quantum 
electrodynamics effects play comparibly impor- 
tant roles. 



II. RADIATION FRICTION EFFECTS 
ON CHARGED PARTICLE MOTION 

In order to describe the relativistic electron 
dynamics in the EM field we shall use the 



equations of electron motion with the radia- 
tion friction force in the Landau-Lifshitz form 
[20[ with a form-factor taking into account the 
QED weakening of the radiation friction. The 
EM wave is modeled by a rotating electric field, 
which as noted above corresponds to the trans- 
formation to the boosted frame of reference 
moving with the group velocity of the wave. In 
this frame of reference the electron equations of 
motion can be written as 



dr 


^rad ^e\Xe) 


7e 


a(q-a) + q {je^f - (q 





7e 



d^ 



if]}, (18) 
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ql is the electron 



where 7 = ^/l 

Lorentz factor. The form factor Ge(Xe) with 

Xe defined by Eq. ([8| describes the radiation 

friction reduction due to quantum effects. In 

3D notation the parameter Xe given by Eq. Q 

reads 




: out, 



vxB -(v-E)^ (19) 



(20) 



Xe = 

We introduce the normalized variables, 

eE 



P 

rripC 



The dimensionless parameter defined by Eq. (pi) 
determines the role of the radiation friction. 

As well known, QED effects weaken the ra- 
diation friction because in the quantum regime 
the charged particle emits less radiation. In the 
QED regime the recoil due to photon emission 
becomes important. According to Refs. [17] 
the total radiated intensity is reduced by a fac- 
tor depending on the quantum parameter Xe- 
In addition, in Eq.dl8| following an approach 
used in Refs. [71 I9( I29[ l3Q] we take into ac- 
count the QED effects by using the form-factor, 
G{Xe)i equal to the ratio of the full radiation 
intensity, /, to the intensity emitted by a clas- 
sical electron. Id- Using the results of Refs. 
[2T1[3TJ[32] we can write the form-factor G{xe) 
as 

Ge(Xe) = 



12 + IbXeX' 



,3/2 



2^3 



uxix 



■^\x)xdx, (21) 



4{l^XeX^/') 



where <l>(x) is the Airy function [33 . Here, we 
neglect the effects of the discret nature of the 
photon emission in the QED (see [2811301 134]). 
In the hmit Xe <^ 1 the form- factor G(xe) 
tends to unity as 



GeiXe) = 1 



55^3 
^6" 



-Xe 



•48x^ 



G 

1.0 



Xp 
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For Xe ^ 1 it tends to zero as 
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327r 
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(23) 



However expression (21) and asymptotical 
dependences (22) and (23) are not convenient 



for implementing them in computer codes. For 
the sake of calculation simplicity the approxi- 
mation 

GsKsiXe) ^ 
1 

[1 + 4.8(1 + Xe)ln(l + l.Txe) + 2.44x2]'/' * 

(24) 
is usually used [29] . This expression can be fur- 
ther simplified. The function 



GRiXe) 



1 



(l + 8.93xe + 2.41x 



2^,2/3' 



(25) 



which we shall use below, within the interval 
< Xe < 10 has an accuracy of approximation 
better than 1%. In Fig. |3]we plot the functions 
Ge(Xe) and Gnixe) given by Eqs. ^ and (|25]), 



FIG. 3. Functions Ge (lower curve) and Gr (upper 
curve) given by Eqs. (21) and (25), respectively, 
versus Xe- 



respectively. We see that that their difference 
is negligebly small. 

We consider a rotating electric field given by 



a = — a (e2 cos r + 63 sin r) , 



(26) 



where 62 and 63 are unit vectors along the co- 
ordinate axis in the plane perpendicular to the 
direction of the EM wave propagation. We rep- 
resent the electron momentum as (see also Ref. 

[35J 




1 

COST sinr 
sin r — cos r 




(27) 



In the case of a rotating electric field Eq. ([19 ) 
yields for the QED parameter 



qI 



il 



as 



(28) 



with 7e 



1 



-Qi 



Ql 



M\' 



Substituting expressions (27) into the equa- 
tions of the electron motion (18), we obtain 



dr 



Q± 



£radGe{Xe)ci q\\ 



il 



(29) 



dq± 
dr 
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adGeiXe) 



leCl- 



a — 1 

7e 



il) 



(30) 



Here and below we do not consider ultrarela- 
tivistic electron beam interaction with the laser 
radiation. We assume that the longitudinal 
component of the electron momentum, qi^ is 
much less than th e transverse momentum com- 
ponent \/<72 + ^3 ' which corresponds to the case 
of the laser pulse interacting with a near-critical 
density plasma, when n ^ rueuo'^ /Ane^ and 

q{ ~ 1, the change of the momentum compo- 
nent is negligible provided that the laser pulse 
duration is small enough. 



= Cl— -eradGe{Xe){ciq^^Cl^q\) , (31) 

7e 



Multiplying equation ( 29 ) by g| | /7e and equa- 
tion ([30| by q^/^e and taFmg the sum we find 

dje_ 

dr 

which shows how the electron acquires energy 
from the electromagnetic wave and loses it due 
to radiation friction. 

Typical solutions of this system of equations 
are presented in Fig. |4J where we show de- 
pendences of the q2 and ^3 components of the 
electron momentum with respect to time time 
(l.h.s column. Fig. |4k,c), and the q\\ and q^ 
components of the electron momentum (central 
coloumn. Fig. U3,d). As we see, for a <C £~^j^ 
the electron osculations in the {q\\^qi_) plane de- 
cay slowly while for the laser pulse amplitude 
values equal to or above e^^^ , the electron os- 
cillations in the rotating coordinate system de- 
cay during a time of the order of or less than 
the wave period. The dependence on time of 
the QED parameters Xe(^) and G{t) is plotted 

-1/3 



in Fig. Wfe,f. We see that for a = 0.25^^ 
i.e. at the laser intensity of the order of ~ 10^^ 
W/cm^ for 1/i m wavelength laser radiation the 
quantum reducing of the radiation friction force 
is negligibly small. Asymptotically at t ^ 00 



have 



Xe 



0.1 and G, 



1/3 



0.95. When 



-1/3 



a = O.TSe^^^^ the quantum correction of the 

the radiation friction force becomes more sig- 

1/3 
nificant with Xe ^ 0-3 and Ge ~ 0.75 in the 

limit t ^ 00. 

Using the above obtained results we can find 

border lines between the domain where the 



dominant radiation friction regime takes place 
(this is the domain II in Fig. |3| and the do- 
main IV in Fig. [3J where QED effects must be 
taken into account. They are given by different 
dependencies of the EM radiation intensity on 
the wave frequency. As was shown in Ref. [35], 
in the limit a > arad the transverse component 
of the electron momentum q± ^ {(i£rad)~^^'^ is 
substantially less than the component parallel 
to the electric field ^n = {a/erad)^^"^- Quantum 
effects become important, when the value of the 
QED parameter (28), which becomes equal to 
Xe ^ Q\\/(^s^ approaches unity. This yields in 
the limit cj <C cji for the EM wave intensity 



^R-Q 



m^c e 
97rn4 
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24 
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(32) 



The border line between domains III and IV 
for CO :^ uj I in Fig. |3] corresponds to the 
limit Xe ^ 1, when the asymptotic expres- 
sion for the form factor Ge(X2) is given by Eq. 
(23). From the equations of the electron mo- 
tion (l29|) and (l30|) we obtain, that the electron 



normalized enery depends on the EM field am- 
plitude as 7e ~ {(^^ /^rad^sy^^' "^^^^ yields 



for the EM wave amplitude at the border line 
aq-R ^ 96as^e;^^ ^ 324(A^/rf A). Using this 
expression for ag-R we obtain a dependence of 
Iq-r{uo) in the form 

Iq-r ^ 



87- 



X 10^ 



UJi 



W 



n2- 



(33) 



In Fig. [5] we plot the border curves 
Ir, Ir-q, Iq and Iq-r versus the EM wave fre- 
quency. These curves subdivide the plane I^uo 
to four domains. In the first domain, (I) neither 
radiation friction nor QED effects are impor- 
tant for the relativistic interaction of an electron 
with the EM field. In the second domain, (II) 
the electron - EM wave interaction is dominated 
by the radiation friction force effects. In the 
third domain, (III) the QED effects come into 
play with insignificant radiation friction force 
effects. In the high intensity limit, (IV) both 
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FIG. 4. Dependence of q2 and q^ (a,b), ^n and q± (c,d) and Ge and Xe (e,f) on time for e 
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as = 4.1 X 10^ (a,c) a = 0.25s;;j^ (d,e) a = 0.75^;^^ 



-1/3 



the QED and radiation friction force effects de- 
termine tlie radiating cliarged particle dynamics 
in the EM wave. The dashed hnes show the de- 
pendences Ir{uj) at < o; < cji and Iq{uj) at 

UJ > UJl. 

The quantum nature of the photon emission 
process results in the recoil effect which causes 
the electron trajectory broadening (see review 
article [36] and literature cited in). Following 
this paper we can write equation describing de- 
pendence on time of the mean square of the elec- 
tron trajectory deflection 



dAr'^ 55 VecXc ^ 

-7e- 



dt 



48^3 A 



(34) 



A condition of the electron deflection during the 
laser wave period 27r /uj being of the order of the 
laser wavelength, Ar ^ A, shows that the quan- 
tum diffusion becomes important at the electron 



quiver energy corresponding to 

1/5 

/ 24V3 mtc'' ] 
Id 



557r huj'^he^ I 



\ 1/5 

a^ \uj) 



2/5 



(35) 

i.e. above approximately 25 GeV. Correspond- 
ing intensity is 



Id = 7.2 X 10^6 (^) 



Wi\4/5 W 



OJ / cm 



2' 



(36) 



III. INTEGRAL SCATTERING CROSS 
SECTION 



According to Eq. (31), the energy flux 



reemitted by the electron is equal to 

e (v • E) ^ 

eradGe{Xe)meC^UJJe (<^<7± + a'^ql) • (37) 





llg^R^O)' / 


IV 


/ ^^ ^ 


^R-0 


/ .y^^ 




/ //^/gGCO) 


" >^ 


I 


^^.ccc,'" 






^ 



For Xe <^ 1, i-e. 



^ 



FIG. 5. Curves Ir{uj) and Iq{oo) subdivide the 
plane I^uj to four domains. I) Relativistic inter- 
action of electron with the EM field when neither 
radiation friction nor QED effects are important. 
II) Electron - EM wave interaction dominated by 
the radiation friction force effects. Ill) QED effects 
come into play with insignificant radiation friction 
force effects. IV) Both the QED and radiation fric- 
tion force effects determine the radiating charged 
particle dynamics in the EM wave. 



This expression is a condition of the balance be- 
tween the rate of energy acquired by the elec- 
tron from the EM field and the rate of the radia- 
tion losses. The integral scattering cross section 
by definition [20] equals the ratio of the reemit- 
ted energy fiux to the Foynting vector ampli- 
tude, cE'^/Att. This yields an expression for the 
cross section dependence on the electron mo- 
mentum 



(T = aTGe{Xe) (<?! + — ) • 



(38) 



Here ar is the Thomson scattering cross sec- 
tion, ar = Sirrl/S = 6.65 x IQ-^^cm^. Expres- 
sion ( 38 ) can also be approximated by 

7e'Ge(Xe) 



(JT 



1 



"rad 



GliXeH 



(39) 



In the range of the wave amplitudes of 1 <C 
a <C {SradGe)~^^'^, wc havc q±_ > ^11 with q^ ^ 
a. The integral scattering cross section grows 
as 



(7t (l 



1 < eJi^ < a < SradCll, 



"rad 



it reaches a maximum of <j 



(41) 



a = l.le 



-1/3 
rad 



0.53cTr£-'i' 



at 



These dependences correspond 



to the domains (I) and (III) in Fig. pi where 
the radiation friction effects are weak. Then for 
a ^ s~^^ the scattering cross section decreases 
according to 



(JT /cis. 



rad-) 



(42) 



as seen in Fig. [6] as a consequence of the fact 
that the maximal power reemited by the elec- 
tron cannot exceed eEc. This gives a constraint 
on the scattering cross section: Aiie/E. In the 



limit £ 



-1/3 
rad 



< a < 



^radO^S 



the electron mo- 



mentum components are q± ^ {ci£rad)~^ and 
(7|| ^ {a/srady^^ with q± <. q\\ [33, which cor- 
responds to domain (II) in Fig. [m where the 
radiation friction is dominant. 

Dependences of the scattering cross section 
normalized by ctt and the electron energy 7e 
on the laser pulse amplitude a are shown in 
Fig. [6) The parameters Srad = 1-75 x 10~^ 
and as = 2.8 x 10^ (curves 2) correspond to the 
EM frequency equal to uji given by Eq. (15). 
For uj = cJi/12.5 (curves 1) the EM wave length 
is about A = 10/im, which is typical for the 
CO2 laser wavelength range. The case 3 with 
UJ = 12.5co'i corresponds to the parameters in 
the domain (IV) in Fig. [6] with the radiation 
friction lowered by quantum effects. 



IV. COMPUTER SIMULATIONS OF 

THE RADIATION FRICTION AND QED 

EFFECTS ON THE LASER PULSE 

INTERACTION WITH 

INHOMOGENEOUS PLASMA 

The electron quivering in the laser field emits 
photons whose energy is proportional to the 
cube of the electron Lorentz factor: hw^ ~ 

0.3fia;7g. For a <C s~^^ a typical value of the 
(40) photon frequency is proportional to a^. In the 




FIG. 6. Dependences of Ig (ct/ctt) and of lg(7e) on 
lg(a). 1) cj = a;i/12.5, 2) cj = cji, 3) cj = 12.5a;i. 
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limit of high laser intensity a ^ e^^^ the fre- 
quency scales as uj^ = ujiajsradf^^ • For multi- 
petawatt laser radiation the emitted photon en- 
ergy is in the gamma ray energy range. 

The gamma-ray pulse energy, duration and 
divergence are determined by the laser pulse 
amplitude and by the plasma target density 
scale length. We analyse the density scale 
length effect on the parameters of the emited 
gamma- flash. Fig. [7| illustrates the concept of 
the high power gamma-ray flash generation in 
the laser-matter interaction. 

During interaction of super-high-power laser 
light with matter the laser pulse is subject to 
various instabilities. Among them the most im- 
portant is the relativistic self-focusing resulting 
in the laser pulse modulation and channeling. 
It leads to the increase of the laser pulse ampli- 
tude and to the decrease of the electron density 
in the interaction region, which change the laser 
energy depletion length and the parameters of 
the gamma-rays emitted. Thorough studying 
of these effects and of the effects of the plasma 
inhomogeneity require computer simulations. 

We performed studies of the laser pulse inter- 
action with high density targets using a two- 



multi-petawatt 
laser beam 





FIG. 7. The concept of high power gamma-ray flash 
generation in the laser-matter interaction. 



dimensional (2D) particle-in-cell (PIC) code 
where the radiation friction force has been in- 
corporated in the Landau-Lifshitz form as has 
also been done in Refs. [37]. In addition, the 
QED radiation friction weakening is taken into 
account by multiplying the radiation force by a 
form- factor Ge(Xe)) given by Eq. (25). We note 
that results of detailed computer simulations of 
the high power gamma ray flash emission by a 
multi-petawatt pulse laser have been presented 
in Ref. [8 , where the QED contribution was 
assumed to be negligebly small. In that paper 
the plasma target inhomogeneity has been sug- 
gested for optimization of the EM wave energy 
convertion to the gamma ray energy flash. 

In the simulations, the laser pulse has the 
power T^ias is equal to lOOPW with the pulse 
duration of Tlas = 30 fs. 

In the tailored plasma target the density 
changes exponentially, n oc exp(x/L), with the 
plasma inhomogeneity scale length, L, from 
0.1 ncr to 350 ncr in the interval ~ 20 /im, and 
then becomes constant having a thickness of 
10 /im. 

The simulation box has a width equal to 
80 /im and a length varying from 50 /im to 
210 /im. The mesh has a spatial resolution of 
Ax = A?/ varying from 1/40 /im to 1/200 /im 
with a temporal resolution of At = 0.0025 fs. 
The plasma is comprised of electrons and ions 
with a mass to charge ratio corresponding to 
AjZ = 2, corresponding to fully ionized Car- 
bon. The average number of quasiparticles is 



10 



about 4 X 10^. 

Simulation results for the parameters of in- 
terest are shown in Figs. [7| and [S] The 100 
PW laser pulse with the normalized ampli- 
tude a = 474 interacts with the plasma tar- 
get, whose density inhomogeneity is character- 
ized by a scale length equal to L = 2.5 jivn. 

Fig. [sja) shows the dependence on time of the 
electron fg^ ion ^i and total energy £^tot when 
the radiation friction force and the QED effects 
are neglected. We see that the laser pulse en- 
ergy at first is converted to the electron energy 
and then to the ion energy. In Fig. Islb) we plot 
the same dependences of the electron, ion and 
total energy and the energy of emitted gamma- 
rays in time for the case with the radiation fric- 
tion force taken into account but without QED 
effects. In this case a major part of the laser 
pulse energy is transferred to the gamma-ray 
flash. Fig. Is] d) presents the gamma ray pulse 
power V^ and energy £^ versus time in this case. 
The maximal gamma ray pulse power is equal 
to 45 PW. If the radiation friction force and 
the QED effects are incorporated into the model 
the conversion of the laser pulse energy to the 
gamma ray radiation becomes lower as it is seen 
in Figs. ItIc) and e). The maximum gamma ray 
pulse power becomes equal to 40 PW. 

Fig. [9] illustrates the hole boring by the laser 
pulse in the tailored plasma targets when nei- 
ther the radiation friction force nor QED effects 
are taken into account (seelola), where the elec- 
tron density distribution in the (x, y) plane is 
shown at t = 140fs and d) with the laser electric 
field distribution presented for the same time), 
when the radiation friction force is present but 
QED effects are neglected ([9] b) and e)), and 
when the radiation friction force and QED ef- 
fects are incorporated into the description ([9] 
c) and f)). The radiation friction makes the 
hole walls less modulated and the laser pulse 
less filamented compared to the case when the 
radiation force effects are neglected. The QED 
corrections lower the radiation force resulting in 
a slightly more modulated laser pulse and hole 
wall. 

The efficiency of the laser energy conversion 
to the energy of the gamma-ray flash, >^e//7 de- 



pends on the laser pulse power. In Fig. 10 
compare the coversion efficiency with and with- 
out the quantum correction of the radiation fric- 
tion for 30 fs laser pulses with the power in the 
range from 0.1 PW to 100 PW. As we see in the 
relatively low power limit the quantum effects 
are negligebly small. At the laser power equal 
to 10 PW the conversion efficiency is 0.25 and 
0.3, and for the 100 PW laser pulse we have 0.55 
and 0.65, respectively. 



V. SPECTRUM OF THE RADIATION 

EMITTED BY AN ENSEMBLE OF 

ULTRARELATIVISTIC ELECTRONS 

The frequency spectrum of a relativistic elec- 
tron rotating along the trajectory with the ra- 
dius of curvature, i?, is given by the expression. 



duj^ 



V3( 



£ 



27r c rripC^ 



00 
u{uj^,£) I K^/s{x)dx. (43) 



i{oj^,S) 



which is well known in the classical electrody- 
namics ^20j. Here Kjy{z) is the modified Bessel 
function [33]. The electron energy is assumed 
to be £^ ::^ rrieC^, rrie and c are the electron 
mass and the speed of light in vacuum, respec- 
tively. The function of the electron energy, 
u{ujj,£) = uj^/ujc^ is the ratio of the ratio of 
the emitted photon frequency u^ to the critical 
frequency equal to 



3c 
2R 



£ 



(44) 



In the limit u{uj^,£) <C 1 expression (43) 
yields dl/duj^ ^ [e^ /c){uj^R/cY^^ , while for 
u{uo^^£) ::^ 1 we have 



dl 

duj^ 



37re2 £ 



2 c rrieC^ 



u{uj^,£fe-<''-^^\ (45) 
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FIG. 8. Computer simulation of a 100 PW laser pulse interaction with a tailored target, a) Time evolution 
of the electron ^e, ion 8i and total energy Stot with no radiation friction force and QED effects, b) 8e, ion 
Si, gamma-ray S^ and total energy £^tot with the radiation friction force but without QED effects, c) 8e, 
ion Si, gamma-ray S^ and total energy Stot with the radiation friction force and with the QED effects, d) 
The gamma ray pulse power V^ and energy S^ vs time with the radiation friction force but without QED 
effects, e) V-y and S^ vs time with the radiation friction force and QED effects. Time is measured in fs. 
Units of the energy and power are J and PW, respectively. 
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FIG. 9. Computer simulation of the hole boring by a 100 PW laser in the tailored target. a,b,c) The 
electron density distribution and d,e,f) and the laser electric field in the {x,y) plane at t = 140fs. For a) 
and d) no radiation friction force and QED effects are taken into account, b) and e) show the hole boring, 
when the radiation friction force is present but without the QED effects, e) and f) correspond to the case 
of the nonvanishing radiation friction force and QED effects. 
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FIG. 10. Comparison of the coversion efficiency of 
the laser energy to the energy of the gamma-rays 
for a 30 fs laser pulse interacting with the tailored 
target when the quantum correction of the radiation 
friction is taken into account (dots) and when it is 
neglected (solid curve). 



tra in the interval 8 < £rn can be approximated 
by a power- law dependence. For the energy 
S < £rn the spectrum has exponential form. 
Under the conditions of the simulations pre- 
sented in Ref. [8 the energy, £rni is of the order 
of 50 MeV. The electron energy distribution can 
be described by the dependence 



dN{8) 
dS 



K£-^exp 






(46) 



For the parameters corresponding to the results 
presented in [8 the power index equals k, ~ 0.8. 
The averaged spectrum of the emitted pho- 
tons is given by the integral 



J{uj^) 



dI^{uj^,£)dN{£) 



duj^ 



d£ 



d£, (47) 



Electrons rotating along the circles with 
the same radius emit electromagnetic radia- 
tion with the same frequency spectrum indepen- 
dently of the particular radiation mechanism. 

In the case of synchrotron radiation when the 
ultrarelativistic electron with the energy £ ro- 
tates in a homogeneous magnetic field the Lar- 
mor radius in the limit £ ^ rrieC^ is equal to 
^B = £/eB. The characteristic photon fre- 
quency of the synchrotron radiation according 
to Eq. (|44| is given by ujc = {3/2)ujB{£/meC^)'^, 



where ujb = eB/rrieC is the Larmor frequency. 

The characteristic energy of a photon emit- 
ted via nonlinear Thomson scattering of the 
circularly polarized electromagnetic wave scales 
with the electron quiver energy as Uc = 
{3/2)ujo{£ /rrieC^)^ ^ where uq is the laser fre- 
quency. The energy of the electron quivering in 
plasma under the action of an electromagnetic 
wave with an amplitude of a = eE/rrieUJC :^ 1 
is of the order of £ = arrieC^. For a laser fre- 
quency of the order of lO^^s"^ the emitted pho- 
ton energy is in the gamma-ray range if a > 10^ 
which corresponds to a laser intensity higher 
than 10^2 W/cm^. 

Computer simulations of a petawatt power 
short laser pulse interaction with a near-critical 
plasma [8 show that the electron energy spec- 



where the functions dLJduo^ and dN{£)/d£ are 



determined by Eqs. (43) and M6^ 



We analyze the spectrum of the radiation 
emitted in the process of nonlinear Thomson 
scattering of a circularly polarized electromag- 
netic wave, when the function u{uj^,£) is 



uc{uj^,£) 



2 Uj f rrieC 

3uJo 



£ 



(48) 



Substituting the functions (43) and (46) to 



(47) and performing integration over £ we ob- 



tain the averaged spectrum of the emitted pho- 
tons 



Jc{ujj) =Qc^-f^ Fc{em,i^), (49) 



where 



V3Ke^ /J_ 



27r c \3ujQ 



{rrieC ) 



2x(l-/.) 



£m f3uJo\ ■ 



and 



rriec'^ \2ujj J 



(50) 



(51) 
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OO 

-(i+«) / 



K^/:i{x)V ( -1 



dx. 



(52) 
Here r(a, z) is the incomplete gamma function 
[33 . The function Fc (e^, n) can be expressed 
in terms of hypergeometric functions. 

In the low frequency range, uj^ <C 
(3/2)co'o(^m/^^eC^)^7 usiug the asymptotic rep- 
resentation of the incomplete gamma function 
at x~^ez} -^ 0, 



Considering the limit e^n ^ 0, which cor- 
responds to the high frequency range, u^ :^ 
(3/2)co'o(fm/^^eC^)^7 wc obtaiu 

Fc{em,f^) ~ 



TT _^ / , 1 + 2^ 

-e^ / dxx 6 exp 





A/3, 



(X exp 



33/4, 



.3/4 



(57) 



r(-l-^) + e-(^+^)f^-e-^ + ...,(53) 
we obtain 



Fc (00, n) = :^ r I — :r- I F 



1 + /^ 



6 



6 



(54) 

where T{x) is the gamma function [3^. This 
expression is valid provided tv > 1. We see that 
the spectrum has the power form. 



Jc(cJ^) ^ UJ^^ . 



(55) 



If the power index is greater than 2, /^ > 2, 
the dependence Jc{^j) is monotonically de- 
creasing as uj^/ujQ -^ 00. For 1 < k, < 2 the 
function Jc{oj^) grows with uj^. In this case as 
well, if the power index is less than unity, hz < 1, 
we should take into account that the electron 
energy spectrum is truncated at the high en- 
ergy end. 

Analyzing the asymptotic behaviour of the 
function Fc (e^, k,) given by Eq. (52) as e^ ^ 



00 for /^ < 1 we find the low frequency scaling 



V3e^ 



(l-K) 



r(i 



1/3 



^'sj 



i. e., the spectrum decreases exponentially with 
photon frequency. 



Jc{^-f 



exp 



51a;of^ 



(58) 



As a result we see that the frequency spectrum 
has a maximum at 



14.2 X ujQ 



f 



(59) 



For the parameters of the simulations presented 
in [8 the characteristic energy £rn is of the order 
of 38 MeV. This yields for the photon energy, 
^7,m = ^-i,mi corresponding to the spectrum 
maximum e^ rn ^ ^ MeV. 



VI. CONCLUSION 

With Kilo- Joule lasers high field science will 
enter novel regimes of electromagnetic radiation 
interaction with matter when radiation friction 
force effects result in the high efficiency conver- 
sion of the energy of the EM wave into the en- 
ergy of hard EM radiation in the form of ultra 
short high power gamma ray ff ashes. The en- 
ergy spectrum of the gamma-ray ffash emitted 
by the relativistic electrons has a typical form 
with a maximum which dependence on the pa- 
(56) rameters is given by Eq. (59). 
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